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ABSTRACT 



From a sample of 12 well-observed Type la supernovae, we find clear evidence of correlations between early phase 
(U — B), (V - R), and (V - I) colors and the velocity shifts of the [Fe II] A7155 and [Ni II] A7378 nebular lines 
measured from late-phase spectra. As these lines are thought to trace the ashes of the initial deflagration process, 
our findings provide additional support to the new paradigm of off-center explosions in Type la supernovae, and we 
interpret these correlations as viewing angle effects in the observed colors. We also show that the nebular velocity 
shifts are related to the strength and width of the Ca II H&lK and IR-triplet lines near-maximum light. The evidence 
presented here implies that the viewing angle must be taken into account when deriving extinction values and distances 
in future cosmological studies based on Type la supernovae. 

Key words, supernovae: general — distance scale 



1. INTRODUCTION 

Type la Supernovae (SNe) play an important role in mod- 
ern astrophysics thanks to their great power in extragalactic 
distance determinations , which permitted the discovery of 
the c osmic acceleration (|Riess et al.ll 19981 iPerlmutter et al.l 
Il999f) . Line-of-sight extinction has to be taken into account 
to correct the distance estimate of the SN host galaxy. 
Several empirical methods that make use of SN colors have 
been developed for this purpose and very different extinc- 
tion laws from the Galactic law need to be invoked to re- 
duce the dispersi on in the calibration of SNe la in cos- 
mological studies (jHicken et al.l [20091 iFolatelli et alJl20ToL 
iWang et alJl2009bh . 

Theoretical and observational studies have shown ev- 
idence of an off-c enter deflagration as the initial burnin g 
proce ss in SNe la (iKasen et al.|[20ul iMaeda et alj|2010dlal 



1201 If) . IMaeda et al.l ( 2010aD show a relation between the 
velocity gradient of the Si II lines observed near maximum 
light and the nebular velocity shifts measured from the 
[Fe II] A7155 and [Ni II] A7378 lines, which are thought 
to trace the region where the in itial deflagration occ urs. 
The off-center explosion model of IMaeda et al.l (|2010al) ex- 
plains the large diversi ty in velocity gradie nts near maxi- 
mum light described in iBenetti et al.l (|2005f) , who classifies 
as either high velocity gradient (HVG) or low velocity gra- 
dient (LVG) SNe depending on the velocity gradient of the 
Si II A6355 line. 

IKasen et all (12009ft . IMaeda et ail 1)201 lh . and 
iFolev fe Kasenl ((201 ll ) demonstrate both theoretically 
and observationally that the (B — V) colors and luminosity 



depend on the asymmetry of the explosion and the viewing 
angle. Here we present additional evidence that favors this 
new paradigm of asymmetric explosions in the form of 
strong correlations of (V — I), (V — R), and (U — B) colors 
at very early phases with the nebular velocity shifts at late 
epochs. These correlations could be useful in improving 
SN color calibrations in future cosmological studies. 

2. SAMPLE AND METHODOLOGY 

Our sample consists of 12 nearby SNe with a well-sampled 
light curve, photometry from at least six days before maxi- 
mum light and late-time spectroscopic observations, i.e. at 
least 150 days after maximum, to measure nebular velocity 
shifts. In Table [TJ we summarize the decline rate (Amf 5 ) 
and line-of-sight reddening values of our Type la SNe sam- 
ple. 

We estimated the epoch of maximum light for each fil- 
ter, Amf 5 , and the SN colors at different epochs using 
fifth order polynomial interpolated light curves. We used K- 
corrected photometry for SN 2004eo and S-corrected pho- 
tometry whenever available in the literature. Given that 
we are interested in assessing the relations between op- 
tical colors and nebular velocity shifts around maximum 
light, we estimated host gal axy reddening using the Lira 
relation ()Phillips et al.l [l999h . which is a method indepen- 
dent of SN colors at maximum . To estimate Galact ic ex- 
tinction, we used the maps of ISchleeel et all (11998T) an d 
the extinction law (Ry = 3.1) of ICardelli et al.l ()l989l ). 
which is also used for host galaxy reddening. The redden- 
ing estimates are shown in Table [T] and the corresponding 
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Tabic 1: Sample of Type la Supernovae. 



SN 


Amf 5 


E(B-V) G al 


E(B-V) host 


Ref. 




1990N 


1.04 ± 0.14 


0.026 


0.163 ± 0.066 


L98 




1994D 


1.40 ± 0.15 


0.022 


-0.028 ± 0.046 


P96 




1998aq 


1.06 ± 0.12 


0.014 


0.060 ± 0.025 


R05 




1998bu 


1.04 ± 0.19 


0.025 


0.401 ± 0.017 


J99 




2001el 


1.16 ± 0.11 


0.014 


0.318 ± 0.020 


K03 




2002bo 


1.02 ±0.23 


0.025 


0.437 ± 0.051 


B04 




2002dj 


1.14 ±0.17 


0.096 


0.149 ±0.050 


P08 




2002er 


1.25 ±0.15 


0.157 


0.189 ±0.057 


P04 




2003du 


1.00 ±0.11 


0.010 


0.032 ± 0.020 


S07 




2004eo 


1.31 ±0.16 


0.109 


0.048 ± 0.024 


P07 




2005cf 


1.00 ±0.13 


0.097 


0.138 ± 0.012 


W09 




2006dd 


1.12 ±0.18 


0.021 


0.054 ± 0.019 


S10 





Refer ences. B04. |Benetti et al.l 
19991); K03 . | K risciun as et alJ 
19981 ); P96. iPatat et all (I1996D ; P04 



(2004]); J99. Una et al 
20031b L98. iLira et al 
Pignata et all (|200^~ 



P07. iPastorello et all J2007l b P08. [Pignata et alJ 



R05. iRiess et all 1120051) ; S0 7. iStanishev et all J2007r i; S10 



iStritzinger et alJ l|2010l b W09. IWang et al.l (|2009al 1 
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Fig. 1: Reddening-corrected colors of our sample of Type la 
supernovae. 



reddening-corrected light curves are shown in Figure U We 
adopt ed the nebular velocity shifts (Vneb) from lMaeda et all 
(|2011h and studied their relation with SN colors at different 
epochs. 

Finally, we investigate the effect of Ca II lines, which 
are the dominant features in the U and / bands, in the 
colors of SNe la. To quantify the effects of Ca II in the 
SN colors, we model the available spectra acquired four 
da ys priot to maximum using the code syn±± described 
m iTh^ mas et all (|2011h . and then we computed synthetic 
photometry for our models with and without Ca II lines. 




-2000 2000 -2000 2000 

V„, [km s-'] V„ cb [km S -'] 



Fig. 2: (top) (U — B) vs Vncb at two different epochs, r is 
the linear correlation coeficient. (bottom) same as 
above but for (V — I)o- 



shift of the [Fe II] A7155 and [Ni II] A7378 lines, which are 
suggeste d to originate from t he ash es of the deflagration 
process (jMaeda et a.l.ll2010dlal . l2T)Tlh . 

We first note that in Figure □ the (U - B) and (V - 
I) colors show a large diversity before and near maximum 
light. The (B — V) and (V — R) color curves, on the other 
hand, are more homogeneous at all epochs. We show that 
some of this diversity in colors could be due to the viewing 
angle. 

We now introduce zero subscripts for the reddening- 
corrected colors. To quantify possible correlations between 
unreddened colors and Vncb, we use the linear model: 



(color) =a + (3 x (V neb / 1000 km s -1 ). 



(1) 



In Figure[2j we present the evolution of the (U — B)o and 
(V — I)o colors as a function of Vneb, at -4 and days since 
maximum light, as well as the best-fitting linear model from 
eq. (H|). From the top and bottom panels, it is clear that 
there is a linear correlation between [U — B)o and V„eb and 
a linear anticorrelation between (V— /)o and V nc b- Although 
not shown here, these correlations are clearly present from 
-8 to ±8 days since maximum light, going generally from 
steeper to shallower slopes with time, in agreement with 
the scatter seen in Figure [TJ A weaker, but significant an- 
ticorrelation was found between (V — R)o and Vneb, and 
no significant correlation was found between {B — V)o and 
Vicb- Best-fitting parameters for eq. ((TJ) at selected colors 
and epochs are shown in Table [5J All the correlations found 
are significant, with a very low probability of the data being 
drawn from an uncorrelated population in all cases. 



3. OPTICAL COLORS AND THE NEBULAR 
VELOCITY SHIFTS 

We now study the relation between optical colors and 
vie wing angle in the context of off-cente r ignition mod- 
els dKasen et al.ll2009t iMaeda et alJl2010dlbM . We assume 
that the viewing angle is traced by the nebular velocity 



4. EFFECT OF VIEWING ANGLE IN Ca II LINES 

Given that the Ca II HSzK and IR-triplet are the strongest 
lines in the U and / filters, respectively, it is pertinent to 
investigate whether or not the strength of these features 
could be responsible for the (U—B)o and (V—I)q variations 
and their dependences on Vneb, presented above. For this 
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Table 2: Correlations between unreddened colors and Vneb 



color 


epoch 


C.V 


P 


T 


(U- 


B) 


-4 


-0.537 ±0.015 


0.106 ±0.002 


0.76 


(U- 


B) 





-0.476 ± 0.012 


0.066 ±0.011 


0.69 


(U- 


B) 


+4 


-0.320 ± 0.013 


0.079 ± 0.015 


0.71 


{V- 


I)o 


-4 


-0.295 ±0.014 


-0.097 ±0.015 


-0.79 


(V- 


I)o 





-0.476 ± 0.012 


-0.083 ±0.012 


-0.83 


(V- 


I)o 


+4 


-0.596 ±0.013 


-0.076 ±0.014 


-0.76 


(V- 


R)o 





-0.071 ±0.010 


-0.038 ± 0.008 


-0.78 



References. Fitting parameters for the model described with 
eq. fTJ for selected colors and epochs. The last column contains 
the Pearson correlation coefficient, whose values imply low prob- 
abilities of the data being drawn from an uncorrelated popula- 
tion: 0.029, 0.04, 0.032, 0.004, 0.001, 0.007, and 0.003 for each 
row, respectively. 




Fig. 3: 



4000 6000 8000 10000 12000 
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Spectral evolution at -4 days since maximum as a 
function of nebular velocity shift (V nc h)- The ob- 
served spectra in black and syn±± models in red 
(dot-dashed lines). The spectra are redshift and ex- 
tinction corrected (Ry = 3.1). The dashed lines cor- 
respond to the expected rest-frame wavelength of 
HSzK and IR-triplet Ca II lines and the spectra 
in th e figure corres pond to SN 1994D ([Patat et al.l 
I1996IV SN 2002bo (jlBenetti et all 1200411 . SN 20r)2d i 



(jPignata etaL)l2008fl . SN 2002er (iKotak et al.l2005h . 
SN 2003du (iSollerman et al.l 120041) SN 2004eo 
(iPastorello et al.l l2007lh SN 2005cf dBufano et al.l 
12009ft . 



purpose, we use the [Fe II] A7155 and [Ni II] A7378 line 
velocity shifts as a proxy for viewing angle. 

Figure [3] shows the optical spectra of our sample of Type 
la SNe sbout four days before maximum, organized as a 
function of Vneb- Th e spectra in Fi g ure E] are corrected for 
extinction using the ICardelli et al.l (|l989ft Galactic extinc- 
tion law (Ry =3.1) and the reddening estimates given in 
Table \T\ shifted to z = 0. 

A clear trend in the width and strength of the Ca II 
H&cK and IR-triplet lines with nebular velocity shift is ev- 
ident in Figure [3] For negative values of the nebular ve- 
locity shift, (i.e. se en from the side clo ser to the deflagra- 
tion center in the iMaeda et al.ll2010al model), the Ca II 



lines are weaker, whereas these features become stronger 
with increasing values for the nebular velocity shift, (i.e. 
seen from the oppos ite side of the initial deflagration in the 
IMaeda et ai1l2010al model). 

At first glance, it seems that the strength of the Ca II 
HSzK and IR-triplet could explain, at least in part, the 
color variations with V nc b. For example, as Vneb increases 
and the Ca II H&.K line gets stronger, the t/-band flux gets 
relatively weaker, making the SN look redder in (U — B)q, 
in qualitative agreement with the trend seen in the top- 
left panel of Figure [5J Likewise, as Vneb increases and the 
Ca II triplet gets stronger, the J-band flux gets relatively 
weaker, making the SN look bluer in (V — I)o, in qualitative 
agreement with the trend seen in the bottom-left panel of 
Figure [5] 

To assess quantitatively the effect of Ca II lines in the 
colors of SNe la, we model th e spectra from Figure [3] using 
syn±± (Tho mas et al.l [20TTh . In syn±±, one computes a 
spectrum by specifying the location and optical depth for 
a given set of ions. The input parameters for synH — h arc 
the photosphcric velocity, the optical depth, the e-folding 
length of the opacity profile, the maximum and minimum 
cut-off velocity for each ion, and the Boltzmann excitation 
temperature for parameterizing line strengths. The spectral 
models are very good over all wavelengths except blueward 
of 3800 A, where the spectrum is dominated by a superpo- 
sition of many metal lines. 

From the resulting model spectra, we compute s ynthetic 
colors using the U BVRI filter functions given by iBesselll 
(|1990t ). We then repeat the calculation after removing the 
Ca lines from the model spectra. Figure |4] shows the total 
effect of the Ca II lines in the (U — B) and (V — I)o colors, 
as a function of V nc h- We note that this is the maximum 
possible effect because we then suppress all of the lines in 
the spectra. In (U — B)o, the synthetic colors vary by ~0.3 
mag, which could account for part of the observed variation 
of ~0.4 mag shown in Figure [2l We do not see the regular 
trend with V ne b expected if Ca lines alone were responsible 
for the correlations found in this work (eq. [TJ. However 
given the difficulties in modeling the spectra shortward of 
3800 A it is difficult to separate the contribution of Ca II 
lines in the colors. In (V — I)o, we recover the trend seen 
in Figure [2j although the overall variation in the synthetic 
colors is a factor of four smaller than the observed variation 
of ^0.4 mag. We conclude that the strength of the Ca II 
lines is only partially responsible for the observed variations 
in colors with Vieb, and part of the color diversity is likely 
due to differences in the slope of the continuum. 

As expected, our synthetic {B — V)o and (V — i?)o colors 
are weakly affected by the Ca II lines, which fall at the 
edge of the B and R bands. In (B — V)o, the color becomes 
slightly (^0.05) bluer as V no b increases from -2000 to 2000 
km s , which is caused by the emission component of the 
Ca II H&zK P-Cygni profile. In (V— R)q, the effect amounts 
to only 0.01 mag over the entire range of Vneb- 

From our syn±± models, we also find that to properly 
model the Ca II IR-triplet of the LVG SNe in Figure |1 it is 
necessary to add a detached high velocity Ca II component. 
In the case of the HVG SNe (i.e. SN 2002er, SN 2002bo, and 
SN 2002dj), using high optical depths and e-folding lengths 
with a single C a II component can re sult in good line mod- 
els. As noted in lFolev fe Kasenl (|201ll ). as the ejecta velocity 
increases, the width in velocity space of the line-forming re- 
gion increases and the lines become broader. Therefore, in 
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2000 4000 

V „.b l km S ~'l 



Fig. 4: (U — B)o and (V — I)o color differences from mod- 
els of the spectra shown in Figure [3J and the same 
spectra after removing the HSzK and IR-triplet Ca 
II lines. 



HVG SNe individual lines overlap and can be modeled with 
a single broad Ca II component. 

5. DISCUSSION AND CONCLUSIONS 



The theoretical work of iKasen et al.l (|2009t ) has put for- 
ward the notion of intrinsically asymmetric ignition and 
subsequent detonation in Type la SNe, to explain the 
"peak luminosity-decline rate relation" and color diversi- 
ties for a given decline rate. Strong support for this idea 
was provided bv lMaeda et al.l (|2010aD by successfully relat- 
in g the early Type la velocity gradient diversity described 
in iBenetti et al.l ((2005) with the late-time nebular velocity 
shifts of the [Fe II] A7155 and [Ni II] A7378 lines, interpret- 
ing the nebular velocity shifts as an indication of viewing 
angle for an initially off-center deflagration. 

We have presented clear evidence of correlations be- 
tween the late-time nebular velocity shifts, or viewing angle, 
and the (U—B)q, (V—R)o and (V— I)o reddening-corrected 
colors from -8 to +8 days since B-band maximum, provid- 
ing additional support to the new paradigm of off-center 
explosions. In (B — V)a, on the other hand, we find no 
significant correlation with the nebular velocity shifts. The 
strengths of these color dependences on V nc b imply that 
the ejecta are asymmetric and that they are as important 
as the relations between color and Amf 5 . 



Vm 



color depen dence 

on velocity gradient fou nd in IFolev fe K asen (20 111) from 



We also note that the (B„ 



a sample of 121 SNe la (jWang et al.ll2009bD could not be 
investigated in our sample because of the small number of 
SNe. On the other hand, the correlation between extinc- 
ti on and Aw^r corre cted {B — V) colors a nd Vn e b found 
Maeda et al.l (|201lf) using the relations bv lFolatelli et al.l 



(2010) are recovered in our sample. Our estimate of the host 
galaxy extinction is systematicall y highe r by 0.08 mag than 
the values given in iMaeda et al.l (|201lj) . These differences 
are possibly caused by our use of the Lira relation, which is 
based on the late-time (B — V) colors, while t hey employ the 
(B mn x — V max ) pseudo color at maximum (jFolatelli et al.l 
2010). We do not find any relation between the residuals of 
our (B max — Vmax) colors or extinction estimates and the 
IMaeda et al.l (|201 lh estimates with Amf 5 or V ne b- 

The reader may wonder whether our results are an arti- 
fact of our method for inferring extinction corrections. This 
concern can be ruled out given that the (U — B) vs V ne b 



and (V — I)q vs V^b correlations go in opposite directions, 
so that it is impossible to wash out both relations simulta- 
neously by changing the reddening technique. 

We also show that there is a relation between V ne b and 
the shape and strength of the Ca II H&zK and IR triplet 
lines, which are the strongest features in the U and I bands, 
i ndepe ndently of Amf 5 . This agrees with iTanaka et al.l 
who noted that HVG SNe tend to have stronger 



Ca II features. As shown in the previous section, the effect 
of the Ca II lines in (U — B)o and (V — I)o is significant, but 
their strength is only partially responsible for the observed 
color variations and their dependence on y nc b is difficult 
to assess. A reasonable conclusion is that the continuum 
emission is also affected by viewing angle effects. 

Finally, the evidence presented here implies that both 
Vneb and Amf 5 must be taken into account when deriv- 
ing extinction values based on maximum-light colors, which 
is a critical parameter in deriving distances from SNe la. 
Furthermore, as shown bv lFolev fe KasenI (|201 lh . ignoring 
this effect could introduce significant biases in the determi- 
nation of extinction laws (Rv)- A larger sample of nearby 
SNe la is urgently required to improve our calibration of 
these correlations and understand the influence of viewing 
angle on SN colors. 
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